Phase diagram of the Hubbard chain with two atoms per cell 
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We obtain the quantum phase diagram of the Hubbard chain with alternating on-site energy at 
half filling. The model is relevant for the ferroelectric perovskites and organic mixed-stack donor- 
acceptor crystals. For any values of the parameters, the band insulator is separated from the Mott 
insulator by a dimer phase. The boundaries are determined accurately by crossing of excited levels 
with particular discrete symmetries. We show that these crossings coincide with jumps of charge 
and spin Berry phases with a clear geometrical meaning. 
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The transition between a band insulator and a Mott 
insulator in a ctae-dimensional (ID) model for ferroelec- 
tric perovskitead has been a subject of great interest in 
recent yearsaia. Tbp model describing this transition, 
originally proposed!!] for the neutral-ionic transition in 
mixed-stack donor-acceptor organic crystalsjiilll3 is 

H = - l YfA+i*<*> + H - c -) + A^(-l) l n^ 



At a fixed value of A_ jsxact-diagonalization studies o 
rings of up to 12 sitesB'til and Hartree-Fock calculations! 
found evidence of a transition with increasing U from 
a band insulating (BI) ionic phase to a Mott insulating 
(MI) quasi neutral phase, which could be expected on 
general grounds. Furthermore, the transition point was 
characterized!] as a metallic point, with divergent delocak 
ization. On the other hand, a field theoretical approach,tJ 
valid in the weak coupling limit (A, U) <C t, concluded 
that a spontaneously dimerized insulating (SDI) phase 
(also called bond-ordered wave) intervenes between the 
BI and MI phases. However, due to the limitations of this 
technique the precise extension of this phase remained 
unknown. Very recently, density, matrix renormalization 
group (DMRG) investigjationsLrtJ and a quantum Monte 
Carlo (QMC) approachB found contradictory evidence: 
the two DMRG calculations reached opposite conclusions 
regarding the existence of the SDI phase, while in the 
QMC results the SDI-MI phase transition was not ob- 
served. Thus, the existence of all these conflicting results 
calls for a further investigation of this model. 

In this Letter we clarify this controversy and accu- 
rately determine the whole ground-state phase diagram 
of Hamiltonian Eq. (Q). This is accomplished by the 
combined use of the method of topological tcappilions 
(MTT) (jumps in charge and spin Berry phases )Ejii3 and 
the method of crossing excitation levcls(MCEL) based on 
the con£ar-mal field theory with renormalization group 
analysisES'Ll These methods, briefly explained below, 



are somehow complementary in the sense that while the 
geometrical content of the MTT is more clearly displayed 
in the strong coupling limit, the MCEL is based on a 
weak-coupling approach. A nice feature here is that they 
turn out to be equivalent for this problem, so that the re- 
sults obtained are expected to be valid for all parameter 
values. State of the art diagonalization of rings with up 
to 16 sites are performed to determine the phase bound- 
aries with errors estimated in a few percent of t. 

The Berry phases are calculated numerically from the 
ground state |g($|, 3>j.)) of -ff($|, $;) in rings of even 
number of sites L threaded by fluxes $ ff for spin a. The 
Hamiltonian H differs from H in that the hopping term 
has the form — t^2 ia .(c i+la .Ci a e' 1 ^'' ^ L +H.c.) . One can map 
H with periodic boundary conditions (BC) into H with 
twisted BC {c i+La — e l ^ a c; La ) using the canonical trans- 
formation Cj a — e l ^"/ L Zj a . The charge (spin) Berry 
phase 7 C (7,5) is the phase captured by the ground state 
when it is followed adiabatically in the cycle < $ < 2-7T, 
keeping $ T = $| = $ ($f = — $| = $). Discretizing 
the interval < <& < 27r into N + 1 points $ r = 2xr/N 
(r = 0,N), the Berry phases are calculated using:l!3 

JV-2 

7c, a = - Jim Im{ln[T\{g(<§> r ,±$> r )\g($ r+1 ,±<$> r+ i)) 

iV— »-oo -*■ ■*- 



x (s^jv-i, ±$jv-i)|<?(27t, ±2tt))]}, 



(2) 



where | ff (27r, ± 2tt)) = exp[ ? f £ ; j(n, T ± n jX )]| 5 (0, 0)). 
An important property of j c is that if the system is 
modified by some perturbation, the change in the po- 
larization P] J- Pi is proportional to the corresponding 
change in 7 C H Here P a is the contribution of electrons 
with spin a to the polarization of the system. Similarly, 
changes in 7,,. are related to changes in the difference 
P| — P| between the electric polarizabilities for spins up 
and downO APj ± APj, = eA7 CiS /27r[mod(e)]. A more 
crucial property is that in systems with inversion symme- 
try 7 C and 7,,. can only be either or ■k, which has led to 
the idea that 7 = (7c, 7s) can be used-as a topological vec- 
tor to characterize different phases.t2l Such a possibility 



is clear in the strong-coupling limit t — ► 0, where (usu- 
ally) all particles are localized: one can choose a gauge 
in which all scalar products in Eq. (0) except the last 
one are equal to 1, so that 7 C is determined by the sum 
i=j- X)i=o J( n j1 ~\~ n jl)- For example, if there is one parti- 
cle per site (U — ► oo), it gives i-n(L — 1) = wr [mod(27ri)] 
for L even, and then 7 C = 7r. Similarly, for a Neel state 
it is easy to see that 7 S = it, and for a charge density 
wave (CDW) with maximum order parameter (A — > oo) 
7 = (0, 0). These values are consistent with the changes 
in Pj ± Pj originated by the charge transport of all elec- 
trons with a given spin to nearest-neighbor sites, required 
to change the extreme BI state with 7 = (0,0) to the 
Neel state with 7 = (tt, %). For the extreme MI state 
(U — > 00), whicb_is,a spin-density wave (SDW), we also 
have 7 = (it, 7r).EM3 By continuity, one might expect that 
these values of 7 characterize also the BI and MI phases 
in weak coupling. As explained below, this is confirmed 
by an analysis based on the MCEL. This change in topo- 
logical parameters (which is sharp even in finite systems) 
indicates non-trivial changes in Pj ± Pj characteristic of 
a phase transition. 

We JBnd another phase with 7 = (tt, 0). From field 
theoryO we know that this corresponds to the SDI phase 
with order parameter D — J2j&(~ l) J ( c j+i<x c jff + H.c). 
If wc consider the more general Hamiltonian 

H' = H- (t AB - t) ^(cJ+toCto + E.c.)(nta-n i+m ) 2 

ia 

+ V ^2 n i<? Tk+lff', (3) 

iaa' 

we confirm that ±hp_Spi phase of H', well established in 
previous studies,li3£3~CJ is smoothly connected with that 
of H for tAB — ► oo- Furthermore, the model of Eq. (0) 
with A = and (V, tAB — t) > contains essentially the 
same phases as H, and allows a more detailed study of 
the relation between the MTT and MCEL. For V = 0, 
while DMRG results in chains of 40 sites are unable to 
detect the opening of an exponentially small gap|L3 the 
MTT predictions with L up to 12li3 .practically coincide 
with those of field theory for tAB ~ lE§li3 and with exact 
results for t A B — > O.Eil 

The MCEL is based on the fact that in a confor- 
mal field theory (which ultimately describes the low- 
energy physics of ID systems in the charge and spin 
sectors if they are gapless) the exponent v of the 
long-distance power-law decay of correlation functions 
(0(x)0(x + d)) ~ dr v is given in terms of excitation 
energies related to the operator 0(x) in the finite ring. 
A crossing of appropriately chosen excitation energies 
for different operators indicates a change in the char- 
acter of the dominant correlations at large distance (a 
phase transition) .EjO The relevant excitation enemies 
for H' with A = have been studied by Nakamura.tH In 
particular, in the weak coupling limit it is known that 
there is a Gaussian transition from the CDW to the SDI 
in H' (A = 0), with the charge gap vanishing only at 



the transition point .LZTIL3 This transition is determined 
by the crossing of the lowest states with opposite parity 
under inversion and the same total momentum K = Ti/a 
(a is the nearest-neighbor distance), calculated with pe- 
riodic (antiperiodic) BC if L = An {An + 2)0 We show 
that this crossing coincides with the jump in j c . In Fig. 
1 we represent P($) = (#($, <&)\H'\g($, $)) for the sim- 
plest case with L multiple of four. Minimizing P(<I>)witIi 
respect to $ and K leads to $ = n, K — ir/a, and 
the Berry phases are obtained following adiabatically this 
state. Using Cj a = e l ^"/ L Zj a ^ one sees that while the to- 
tal wave vector K of H' remains constant as $ is changed, 
in general K = K + (iV T $ T + ATj$;)/(i a), where N a 
is the number of particles with spin a . The conditions 
leading to the minimum energy ($ = ■k, K = vr/a) cor- 
respond to antiperiodic BC and K = in H' . 
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FIG. 1. Energy per site as a function of flux $| = $| = $ 
for the two lowest lying eigenstates within the subspace of 
K = 7r/o, S pf„0 1 Ip^ameters are L — 4, A = 0, tAB — t — 1, 
V = 2 and three different values of U as indicated. Full 
(dashed) line correspond to states with 7 C = (fy c = n). For 
U — U C (L) 7 C is undefined due to the degeneracy at <E> = 0. 



It is easy to see that the inversion Ic.fi = c_ ■ 



e l ^ a c* L _ - CT is a symmetry of ff'($|, $x) only if both $ CT 
are either or tt (corresponding to periodic or antiperi- 
odic BC). At U = U C (L) and $ = there is a crossing 
of the lowest levels with K = ir/a and total spin S = 0. 
This crossing is possible because the corresponding wave 
functions have opposite parity, and therefore they do not 



mix at $ = 0. For $ — > one can use perturbation the- 
ory in <f>, and for U near U C (L) the state \g(&)) is deter- 
mined by a 2 x 2 matrix involving the above mentioned 
two states for $ = with off-diagonal matrix elements 
linear in $. From the trivial solution to this problem one 
realizes that the product in Eq. (||) for U — * U C (L) — + 
differs from that for U — > U C (L) + + in sign. Hence, 7 C 
jumps at the same place where the transitions occurs ac- 
cording to the MCEL. While for the t-U-V model (H' 
with A = = tAB — t) in weak coupling the transition 
is second order and the charge gap vanishes at the tran- 
sition, for (U,V) ^S> t the transition is first ordet^f 
a fully gapped CDW to a charge gapped SDW.EZrE; 
As a consequence, the MCEL looses its support from the 
conformal invariant (massless) theory. However, in this 
limit the geometrical meaning of the jump in j c is very 
clear, as explained earlier, and justifies the method. 

In the MCEL, the Kosterlitz-Thouless transition, 
which corresponds to the opening of a spin gap, is de- 
tected through the crossing of a singlet even under in- 
version with an odd triplet, with K = and periodic 
BC (for L = 4n)!B At $ = or $ = n, H' ($,-$) 
has SU(2) and inversion symmetries, which are lost for 
other values of $. Therefore, a similar analysis as above 
shows that 7 S jumps at this point JIj If L = An + 2, peri- 
odic and antiperiodic BC are interchanged. When A ^ 
the symmetry under translations in one lattice spacing 
a is lost, K = n/a becomes equivalent to K — 0, and 
the CDW order parameter is different from zero also in 
the SDI and MI phases. The field theory results for H 
show that for XA, U) <t the spin transition retains the 
same featuresjj The charge transition, which for A = 
is described by the sine Gordon model, for A ^ is 
determined by the double sine Gordon model, and the 
universality class changes from Gaussian to Ising. How- 
ever, the transition remains second order and the charge 
gap vanishes at the transition^ Then, at this point and 
sufficiently low energies the charge sector is described by 
a conformal invariant theory, justifying the MCEL. 

In spite of the breaking of traslational symmetry, fortu- 
nately the relevant crossings for H can still be identified 
looking for the ground state energy in subspaces with 
K = 0, total spin projection S z — 0, and a definite par- 
ity under inversion and time reversal. The latter allows 
us to separate states with even and odd S. If the more 
general model H' with A ^ is considered, there are 
some regions of parameters in which the charge transition 
corresponds to a crossing of first excited states within 
the above mentioned subspaces, but we restrict ourselves 
here to the phase diagram of H. For this model, the con- 
nection between the jump in j c and a symmetry switch 
of the gpcwnd state for appropriate BC has been noted 
earlier pun but the relation with the MCEL has not been 
discussed. Moreover, neither results for 7 S nor numerical 
investigation of the SDI-MI transition has been reported 
so far. The calculation of 7 S in H presents technical dif- 
ficulties due to additional crossing of levels (not related 



with phase transitions) which take place for K = 0. We 
have verified numerically that the jumps in 7 C and j s 
correspond to the above mentioned level crossings. 

For given A, we have calculated the critical on-site 
repulsion U c (U s ) at which the charge (spin) transition 
takes place. In addition, for small U and A we have fixed 
U and determined the critical values A c and A s . This 
was done by fitting a quadratic polynomial in 1/L 2 to 
the results for L = 10, 12, 14 and 16, followed by an 
extrapolation to L = oo. This fit works very well for 
A > 0.25 (we set t = 1 as the unit of energy unless 
otherwise stated), and improves with increasing A. The 
difference between U S (L) and U s (L + 2) rapidly decreases 
with L if A is not too small. Instead, for small values of A 
the finite size effects increase and, as a consequence, the 
error in the extrapolation becomes larger. To estimate 
this error we have repeated the fits using L = 8, 10, 12 
and 14; for A = 0.05 this gives a new estimation of U c 
(U s ) that differs from the previous one in 0.12 (0.07). 
For A < 0.05 the relative error in U c and U s becomes 
very large, and we do not present results since they loose 
quantitative validity (except at A = 0, where U c — U s — 
for all L). Instead, for A > 0.25 the estimated error in 
U c , U s is less than 0.06, and less than 0.03 for A > 0.5. 

The resulting phase diagram is presented in Fig. 2. In 
qualitative agreement with field theory results, and for 
all values of A, the transition from the BI phase to the MI 
phase with increasing U takes place in two steps: first, 
a charge transition to the SDI phase at U = U c occurs, 
and then, for U — U s > U c , the spin gap closes. The 
behavior in the strong coupling limit is quite different 
from that of the t — U — V model, for which a first order 
CDW-SDW transition occurs and is_easilv understood 
in terms of perturbation theory (PT)EZIc3 , Ej. Instead, H 
remains non-trivial for t — > as long as U — 2 A ~ t, since 
charge fluctuations are still possible. As a consequence 
of this derealization of charges, 7 C inside the SDI phase 
cannot be calculated analytically just adding the position 
of the charges, as we explained before for A — > oo or 
U — > oo. For t = the SDI phase is absent, and PT in 
t diverges at U = 2A where the BI-MI transition takes 
place. For t < \U - 2A| < U PT is valid, and can be 
used to calculate the energy of the BI and MI phases for 
negative and positive U — 2 A respectively. The MI phase 
in this limit is described by a Heisenberg model with 
exchange J — 2t 2 /(U — 2A). The energies up to second 
order in t are given by E B i = U - 2A + At 2 /(U - 2A) 
and Emi — —J In 2. While the SDI phase cannot be 
described by PT in t, its boundaries are very accurately 
determined by our method for small t. The jumps in 
Berry phases have very little size dependence and show 
that U c ~ 2A + 1.33£and[/ s ~ 2A+1.91tfor t < (A,j7). 
The fact that the SDI phase exists for positive values of 
U — 2A was to be expected from the asymmetry of Ebi 
and Emi under a change of sign of U — 2A. 

The results for (A, 17) » t can be extended qualita- 
tively to A ~ t. The SDI has a nearly constant width 
~ 0.6i, and both U s and U c increase slightly with de- 



creasing A. For A < 1 the critical values U c and U s de- 
crease abruptly, until they reach U c = U s = at A = 0. 
However, in the region < A < 0.25 (0 < U < 2) the 
relative errors in U c , U s become larger with decreasing 
A; in particular, for A < 0.1 our results are not quanti- 
tatively reliable. For (A, U) <C t, the spin transition can 
be estimated integrating out the charge degrees of free- 
dom, assuming that they are described by a free massive 
boson. This leads to a renormalization of the effective 
interaction gi± responsible for the opening of a spin gap 
in the sine Gorden model which describes the spin sector 
at low energies B From the vanishing of the renormal- 
ized gi± one obtains the approximate field theory result 
A{* ~ E g ^U/{8irt), where E g is the gap for A = and 
is known from the Bethe ansatz solution. For U -C t, 
Eg = (8/ir)y/tUexp(— 2irt/U), and the exponential de- 
pendence dominates the behavior of A{*. Due to the 
numerical uncertainties for A < 0.1, we cannot establish 
where this exponential dependence deviates from the ac- 
tual SDI-MI boundary. For small A and any value of 
U, an accurate field theory result for A s might be ob- 
tained using a bosonization aaproach which starts from 
the exact solution for A = 0.E3 
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FIG. 2. Ground-state phase diagram of H at half filling. 
The dashed line corresponds to the field theory result A{* (see 
main text) . The open (full) symbols were obtained keeping U 
(A) constant. 

The SDI order parameter D couples directly with op- 
tical phonons with wave vector K = and, therefore, the 
latter should increase the extension of this phase. In prin- 
ciple, one can include these phonons in the numerical cal- 



culations using the adiabatic approximation. However, 
due to the breaking of inversion symmetry our method 
cannot be used to find the phase boundaries in this case. 
QMC calculations suggest that in the adiabatic approxi- 
mation the whole MI phase disappears and the SDI takes 
its place.El This is not necessarily the case if the dynam- 
ics of the phonons is includedfl We must emphasize that 
in the MI phase both dimer-dimer and spin-spin corre- 
lation functions have the same leading power-law decay 
at large distances. The dominance of spin-spin correla- 
tions dtte-tp-logarithmic corrections characterizes the MI 
phase jlalll'ta This renders it very difficult to determine 
the SDI-MI boundary-by direct numerical evaluation of 
correlation functionsJia 

In summary, we have determined the quantum phase 
diagram of the Hubbard chain with alternating on-site 
energies at half filling using topological transitions. The 
method is justified from geometrical considerations in the 
strong coupling limit (t — > 0) and by field theory argu- 
ments in the weak-coupling (U, A) <C t region. We con- 
firmed the existence of a spontaneously dimerized phase 
and determined its boundaries for the first time. 
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